A probability-conserving dissipative Schrodinger equation 
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Dissipative eflects on a microscopic level are included in the Schrodinger equation. When the 
decay between different local levels as a result of the coupling to a bath, the Schrodinger equation 
no longer conserves energy, but the probability of the states is conserved. The procedure is illustrated 
with several examples that include direct electronic decay and damping of local phonons (vibrational 
levels). This method significantly reduces the calculational effort compared to conventional density 
matrix techniques. 
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Introduction.— A very diverse range of research fields 
such as solar energy conversion, quantum information 
storage, photosynthesis, and light controlled devices de- 
pend on a deeper understanding of the nonequilibrium 
dynamics of the system under consideration. The in- 
creased use of pump-probe type experiments, in partic- 
ular with the coming on line of the next-generation X- 
ray sources, provides increasingly detailed knowledge of 
the transient electronic properties of a wide variety of 
materials [1[. The conventional approach to nonequilib- 
rium problems is the density matrix method [2|- The 
diagonal elements pu{t) of the density matrix describe 
the population or probability of state j as a function of 
time t, whereas the off-diagonal matrix elements pui (t) 
with i' ^ i give the phase or coherence of the system. 
This method is a powerful approach when one wants to 
obtain a detailed understanding of the nonequilibrium 
dynamics. However, the density matrix method is not 
entirely without problems. Within certain approxima- 
tion schemes, the density matrix may become diverging 
[3| . This problem can be avoided by the use of Lindblad 
operators [j|, whose implementation is often far from 
straightforward, even for simple systems. Although this 
can in principle be solved, a more serious disadvantage 
of the density-matrix technique is the difficulty in han- 
dling large systems. In many problems, the size of the 
basis TV needed to describe the local system can easily 
be of the order TV ~ 10^~^. Within a density matrix 
framework, this requires the evaluation of N"^ density 
matrix elements, which is intractable even for small TV 
and impossible for large N . In contrast, if one could 
directly study the nonequilibrium dynamics of the wave- 
function IV'(i)) = X]i=i Ci(i)IV'i): where the ji/'i) form the 
basis of the system under consideration in bra-ket no- 
tation, only N coefficients need to be calculated. The 
density matrices can then be directly obtained through 
p,^{t)=c\{t)c,(t). 

The goal of this Letter is to describe dissipa- 
tive processes with an effective probability-conserving 
Schrodinger equation. 



ih 



d\mi 

dt 



where Hq is the Hamiltonian of the system and D de- 
scribes the effective dissipation. Several phenomeno- 
logical expressions for D have been suggested, such as 
D - (l/n)ln(^/V'*) @ and i:> - {d / dt) In^iP* iP) 0. 
However, these damping terms primarily dealt with the 
quantum behavior of macroscopic variables, such as a 
Ginzburg-Landau type wavefunction or the phase differ- 
ence across a Josephson tunnel junction. Our interest 
lies in incorporating microscopic dissipative mechanisms 
within systems consisting of, for example, an ion (e.g. 
transition metal or rare earth) surrounded by ligands. 
This "local" system is considered part of a larger system 
such as a molecule in solution or a solid. The latter con- 
stitute the effective surroundings that can dissipate en- 
ergy from the local system. First, we describe the nature 
of the dissipative Schrodinger equation in the presence of 
a microscopic dissipation mechanism. Our focus lies on 
ultrafast intersystem couplings and we only consider the 
zero-temperature formalism. We then demonstrate how 
to incorporate electronic and vibronic dissipation. 

Dissipative Schrodinger equation.— The Schrodinger 
equation in the presence of a bath is given by 
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where Hq is the Hamiltonian of the local system and Hb 
the Hamiltonian of the bath plus the interaction of the 
local system with the bath. When writing the coefficient 
Ci for a particular vector in the basis in terms of an ampli- 
tude ai{t) = \ci{t)\ and a phase fi, or Ci{t) = ai(i)e*''''(*\ 
we can write for the change in coefficient due to the pres- 
ence of the bath 
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The latter terms gives the change in phase, which causes 
an embedding of the local system in its surroundings. 
Since the bath has many degrees of freedom, we assume 
that the phase change of the local system due to the bath 
can be neglected based on law of large numbers. In ad- 
dition, due to the complexity of the surroundings, the 
precise nature of this embedding is often very difficult 



and can usually only be taken into account in some effec- 
tive way. We therefore consider only the strength of the 
decay. Now let us assume for the moment that we are 
able to determine an expression for the change in ampli- 
tude Pi = af in a particular basis (e.g in the absence of 
intcrsystem coupling constants) 



dt 



2a,; 



dai 
It 



I{{P^}l 



(4) 



where / is a function of the probabilities Pi . In the follow- 
ing, we demonstrate how to derive differential equations 
for Pi and consider the effects of their incorporation. The 
change in the coefficient due to the bath is then given by 
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This leads to a dissipative term in Eq. ([T} given by 
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The dissipation does not necessarily have to be diagonal. 
After deriving the diagonal dissipation in a particular 
basis, a unitary transformation to a different (more suit- 
able) basis can be made. Note that if Eq. Q conserves 
the total probability, the probability is also conserved in 
the dissipative Schrodinger equation. 

Electronic transitions.— In the following, we give sev- 
eral examples that demonstrate the use of the dissipative 
Schrodinger equation. First, we consider electronic tran- 
sitions, where a decay mechanism can be given that di- 
rectly couples the states of the local system. This decay 
mechanism does not necessarily conserve the number of 
electrons of the local system, i.e. the problem can also be 
classified as an " open quantum system" . A typical ex- 
ample is the Fano problem [3| , where a local state is cou- 
pled to a continuum and the particle number at the local 
state is not conserved. However, alternatively, one can 
also view open quantum problems as closed systems when 
considering the probability of the local states, which is 
is conserved. In the Fano problem, there are two states: 
|1) and |0) with and without a particle on the local site. 
The total probability for the local states is conserved and 
governed by the equations 
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The corresponding dissipative Schrodinger equation gives 

(8) 
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The solution is straightforward and gives, as expected, 
ci — e^'"* and cq — \f\^-e^^ . 

Although the result for the Fano-type problem is 
hardly surprising, the real advantage lies in situations 
where both electronic coupling and decay are present. 
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FIG. 1: (color online) The time dependence of the probabili- 
ties of a four-level system, with hybridization between states 1 
and 2 and a cascade from 2 — >■ 3 ^- 4, see inset. {l|_ffo|2) = 0.3 
eV and {2|ii-o|2) - (l|i/o|l) = 0.5 eV. dPa/di = -IV^zPi, 
dPs/dt = -2r34P3 -I- 2r23-P2, and dP^/dt = 2r34-P3 with 
ftr23 = 20 meV (corresponding to a decay time of 66 fs), 
fir 34 = 10 meV (33 fs). 



Figure [T] shows an example with four levels, where the 
details are given in the caption. The initial state 1 is cou- 
pled via hybridization to level 2, which is higher in energy 
(the parameters are given in the caption). State 2 then 
shows a cascade decay via level 3 into level 4. State 1 and 
2 show the typical quantum oscillations associated with 
hopping between two levels. Some oscillations are still 
observed in state 3 due to its coupling to state 2. The 
probability of finding state 4 steadily increases as that 
of initial state 1 decreases through intermediate states 2 
and 3. 

Phonon decay.— In many dissipative processes, the 
number of particles in the local system is conserved and 
energy is lost through the damping of atomic vibrations. 
This situation occurs, for example, when two states with 
different interactions with the surroundings are coupled. 
Following an intcrsystem crossing, phonon modes (cor- 
responding to, for example, local vibronic oscillations) 
are excited that decay on the femtosecond scale through, 
e.g. intramolecular energy redistributions. Let us first 
consider the decay of the excited phonon states and then 
incorporate this dissipative mechanism into a model for 
intersystem crossing. 

Let us take a Hamiltonian that couples a local vibra- 
tional mode to a bath of phonon modes 



2rPi(i). (7) Hb = nuja'^a + Yyy^kb\bk + Vk{abl + ahl + Yi..c.)], 
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where a^ and h], are the creation operators for an exci- 
tation in the local and continuum modes of energy hw 
and fiwk, respectively. The last term with Vk describes 
the quadratic coupling to the bath. We want to obtain 
an equation for the probabilities P„ of having n excited 
phonons. The fastest damping processes are primarily 
due to intramolecular energy redistributions, which are 
generally almost independent of temperature. We there- 
fore take the limit of zero temperature. The damping 




400 600 800 1000 

afs) 

FIG. 2: (color online) Exact solutions of phonon damping 
for a small number of excited phonons. (a) Damping for 
an initial state with n — 1 phonon. The red lines corre- 
spond to the probability of the effective Ticff = 1 consisting 
of a^|0) (dashed) which couples to the a^a^6J,|0) states (dot- 
ted gives the total probability). The solid line gives the sum 
of both probabilities. The blue lines correspond to the ef- 
fective Tieff = states consisting of the bj^|0) (dashed) and 
the a'bl,by\0) (dotted) states. The solid blue line gives the 
sum of both probabilities. The thick dots correspond to the 
rotating- wave approximation (only the o^lO) and 6j,|0) are in- 
cluded). Note that in the limit F <^ fiu, the results are almost 
identical, (b) Phonon damping for n = 2 in the rotating wave 
approximation. The probabilities for n = 2 (red), 1 (green), 
and (blue) are shown. The dashed lines give the analytical 
results. 



due to Hb can be solved numerically, but only for a very 
limited number of excited levels. Figure [5] shovifs the case 
for n — 1 for a local vibronic mode vifith huj = 30 meV 
and a continuum of bath levels from to Shu: the cou- 
pling Vfc = 6 meV. The contribution of the terms that 
do not conserve the number of phonons, a^fej, and abk, is 
relatively small since these excitations have an energy of 
at least huj. Particle conserving terms can occur at the 
same energy. Furthermore, the non-phonon conserving 
terms act as an effective embedding. The effective boson 
a\l) + ^j, /3fe|2, k), where \n, k) is a state with n excita- 
tions of the local vibronic mode and a boson with index 
k in the continuum, decays to a'\Q,k) + X^fefc' /^fcfc'l-'^'^) 
in an almost identical fashion as |1) decays to |0), see 
Fig. [5] This is only correct in the limit that the broad- 
ening of the local vibronic level due to the continuum 
is less than huj. However, since an intersystem cross- 
ing on a timescale of 100 fs corresponds to an effective 
broadening of 3 meV, this condition is generally satisfied. 
In the remainder, we assume that the ultrafast decay 
can be described by (effective) bosons that couple to the 
continuum in a particle-conserving manner. This is also 
known as the rotating-wave approximation. We want to 
study the equation of motions for a system with n ex- 



cited local vibronic modes A]^|0), where A^ = (a''')"/vnl 
and |0) is the vacuum state. The time dependence of 

Al{t) = e-iff«/?»^1^e*fft/'i is given by 

^^ = -^nu:Ai{t) - ^V^Y. ^^^I W^-iW- (9) 

k 

A similar equation can be obtained for ^{.(i). Using a 

Markov-type approximation, X^fc ^fc^ /n '^^'^^*'^'°^* * ^ ^ 
T5{t-t' 

Aiit) = 



the solution can be obtained 



(10) 



-le 



-inzt 



'n 



fdt'j2Vk^^"^'-''^''biAi_,{t'), 

■Jo ,. 



with the complex frequency z — u — iV. For A'^ = 1, the 
problem is equivalent to a Fano-type decay and we can 
write A\{t) = [ci{t)a'^ + Ejt Cfc(0^l]|0)- The probabili- 



ties are given by pi{t) = |ci(i)|^ 



-2Tt 



and poit) 



J2k \^k{t)\ = 1 — e~ . After obtaining the result for 
A^ = 1, we can derive general expression for arbitrary A 
via iteration. This problem is separable and the solution, 
after a somewhat lengthy but straightforward deriva- 
tion, is given by A^t) = [ci{t)a^ + EkCkit)bl]^\0). 
The probability of finding n local phonon modes at a 
particular time for an initial condition of A^ phonon, 
P„{t) = (A^!/(A - ny.n[)p'^p^-'\ In general, an equa- 
tion can be derived for the probability of n phonons 



dPnjt) 

dt 



-2nTPn{t) + 2{n + l)rP„+i(i), (11) 



which is the microscopic decay equation Q that can be 
incorporated in the dissipative Schrodinger equation. 
We now include the decay in, e.g. the Hamiltonian [3] 



H 



fiLL!a^a+ y [[Ei + \/eihLL!{a + a^)ni] 

i=l,2 

+V{c\c2+clci), (12) 



which describes two states at energy Ei that couple dif- 
ferently to a phonon (e.g. a local vibronic mode) of en- 
ergy huj via the coupling \/eihuj; the states couple with 
each other with a strength V. Since states 1 and 2 have 
different equilibrium positions, a transition between be- 
tween induces an oscillation of the vibronic mode. Using 
the dissipative Schrodinger in Eqs. ^ and (|6]) in combi- 
nation with the change in phonon occupations as given 
by Eq. (fTTj) . we can effectively damp these oscillations. 
However, we need to make a choice in what basis to de- 
fine the decay. In the limit V = 0, the Hamiltonian 
separates into two independent displaced harmonic os- 
cillators. Defining the occupations P„ in the displaced 
bases leads to the physically intuitive situation that, in 
the absence of V, each state relaxes to its equilibrium po- 
sition. Since, only the relative equilibrium positions are 
of importance, we can take £i = and £2 = £■ Account- 
ing for the phonon self-energy Si, the bias is then given 
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FIG. 3: (color online) (a) Occupation Pi of the lowest vibra- 
tional level (n — 0) of state 1 as a function of time. There is 
no bias between states 1 and 2 (A = 0), see inset. The initial 
state is state 1 with n — 0. The dotted blue and solid red 
lines give the result in the absence and presence of damping, 
respectively, (b) The occupations Pi of the lowest vibrational 
levels of states 1 (red) and 2 (blue) in for a bias equal to the 
phonon self-energy (A = e), see inset. The dotted and solid 
lines are in the absence and presence of damping, respectively. 
Note the different timescales for (a) and (b). The parameters 
are e = 0.2 eV, &j =30 meV, V =20 meV, KT = 10 meV. 



hy A ^ El — E2 + e. Figure [3] demonstrates the effects of 
including the damping of the phonon oscillations on the 
occupations. The parameters are given in the caption. 

In Fig. ^a), the time-dependence occupation of the 
lowest vibrational level (n — 0) of state 1 in the absence 
of a bias (A = 0, sec inset in Fig. EJa)) is given. This 
state is also the initial wavefunction. In the absence of 
damping of the local mode, there is an oscillation between 
the n — levels of state 1 and 2. Since the direct coupling 
between n = levels of state 1 and 2 is Ve~'^^^^ is very 
small, the effective coupling between the two levels oc- 



curs through the states where the Franck-Condon factors 
are large. This effectively creates an energy barrier with 
a height of the order of e between the two states. The 
inclusion of damping causes a decoherence of the oscilla- 
tions that damp out jlfll • Figure ^h) shows the results 
in the presence of a bias (A — e), see inset. For this 
value of the bias, the coupling is relatively strong since 
the lowest vibrational level of state 1 is at the maximum 
of the Franck-Condon continuum 111. In the absence of 



dissipation, there is a periodic recurrence of occupation 
in state 1, related to the oscillation period of the local 
phonon. The occupation of the lowest vibrational level 
of state 2 is negligibly small in this case. When damping 
of the local phonon mode is introduced, this recurrence 
is strongly reduced [l2] • In addition, there is a steady in- 
crease in occupation of n = of state 2, which becomes 
close to 100% occupied on the timescale of a few 100 fs. 

Conclusions.— We have demonstrated how dissipa- 
tive effects can be included on a microscopic level in a 
Schrodinger equation by incorporating differential equa- 
tions that describe the change in microscopic occupa- 
tions. The results are intuitive and relatively straight- 
forward to incorporate. If the decay equations are prop- 
erly defined, the Schrodinger equation does not conserve 
energy, although probability is still conserved. The pro- 
cedure was illustrated with examples that included a di- 
rect electronic decay and a damping of local phonon lev- 
els. The described method makes the size of the calcula- 
tional efforts proportional to the size of the basis set N, 
which is a strong reduction to the commonly-used den- 
sity matrix method (proportional to N'^). This makes 
this approach extremely useful in describing relatively 
large systems that undergo ultrafast decay of an excited 
state. Future research should involve the inclusion of the 
effects of temperature. 
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